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^>«^ ' We calculate the electronic structure of nm scale InP islands embedded in Gao.51Ino.49P. The 
\ calculations are done in the envelope approximation and include the effects of strain, piezoelectric 
polarization, and mixing among 6 valence bands. The electrons are confined within the entire island, 
while the holes are confined to strain induced pockets. One pocket forms a ring at the bottom of 
^ ', the island near the substrate interface, while the other is above the island in the GalnP. The two 
■ sets of hole states are decoupled. Polarization dependent dipole matrix elements are calculated for 
' both types of hole states. 
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I. INTRODUCTION 



In recent years great progress has been made in the fabrication and measurement of semiconductor quantum dots 
^ . made by Stranski-Krastonov growth. In this technique, material is deposited epitaxially onto a substrate that is 
■ lattice mismatched to the deposited material. Due to the mismatch, the deposited material spontaneously forms nm 

a scale islands, the size and shape of which are both material and growth-condition dependent. The pyramidal shape 
I ' and presence of strain in the islands makes for a potentially rich electronic structure. Theoretical studies of strained 
'"^ islands have employed various degrees of approximation to the geometry, strain distribution, and electron dynamics, 
^ ■ ranging from single band models of hydrostatically strained islands, to multiband models including realistic shapes 
O \ and strain distributions [|l]-Q. 

I In this paper we consider InP islands bm'ied in Gao.51Ino.49P (GalnP) a material combination that has so far 

received little theoretical attention. The geometry is made as realistic as possible by adopting the shape shown in Fig. 
CO 1, which has been observed using atomic force microscopy and transmission electron microscopy [||. We assume that 
^ different size islands have the same shape as shown in Fig. 1, but are simply rescaled. We parameterize the size of the 
, island by the height in the z direction, h. In most of the calculations we set ft = 15 nm, but we will consider variations 
' as well. Note that the wetting layer is not included in our calculation. The wetting layer does not significantly alter 
ly-j ' the strain distribution in the island because it is pseudomorphically strained. Also, its measured thickness is on the 
order of the grid spacing used, and hence its inclusion in the electronic structure portion of the calculation would be 
' unjustified. The calculation is done in two steps: the strain distribution is determined using continuum elastic theory, 
On ■ and then used in a strain dependent k • p Hamiltonian. 



II. STRAIN 



rH ' The strain is computed by minimizing the free energy for a cubic crystal, 

o ■ 



C xxyyip^ xx^ yy H~ ^xx^zz ~\~ ^yy^zz^~^ 



^ ; 2Cxyxy{(T^y + cr^^ + (TyJ - a{axx + <Jyy + Cr^^), (2.1) 

where (Xy is the strain tensor, the C"s are material depended elastic constants, and a is a parameter used to enforce the 
lattice mismatch between the two materials. The strain is given in terms of the displacement by (Ty = ^{diUj +djUi). 
F is constructed as a function of the UiS on a cubic grid with derivatives replaced by differences, and then F is 
minimized using the conjugate gradient algorithm. The grid consisted of a periodic box of 130 x 130 x 120 sites with 
the island contained within a 65 x 65 x 20 region. The grid is reduced 1/4 by exploiting the symmetries of the system. 
The size of the box is chosen sufficiently large to make the band energies within 1% of the bulk GalnP values at the 
points furthest from the island. Errors due to the finite system should be < 1% at the island. No constraints need to 
be directly imposed on the minimization since far from the island the strain has relaxed to its bulk (= 0) value. 
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The lattice mismatch between the two materials is implemented as follows. The grid is chosen to be commensurate 
with the barrier material so the grid corresponds to some multiple of the barrier material's crystal unit cell. Then for 
unstrained barrier material we would have = 0. But this means that the grid is wrong for the island material since 
it has a different crystal unit cell. Unstrained island material should actually have a nonzero Ui to reflect the fact that 
its crystal unit cell has a different size. If the island and barrier lattice constants are a/ and as respectively, then 
setting a — {Cxxxx + '2Cxxyy){0'i — clb)/0'B in the island causes the minimmn of F to be at a^^ — Sij{aj — aB)/aB, 
which is the desired amount of strain in the grid. Because af^ is completely fictitious, and reflects only the straining 
of the coordinate system, the physical strain is obtained by subtracting cr^^ from the (7^ that minimizes F for the 
heterostructure. 

There is some freedom in the choice of grid differences used to approximate F . Using symmetric differences, 
dxf{x) — > [f{x + h) — f{x — h)]/2a, produces oscillatory solutions that must be smoothed. Oscillatory solutions are 
in fact a common problem with finite difference approximations to equations involving first derivatives. For the case 
at hand this problem is avoided by evaluating F within each unit cube as a function of the u^'s at the eight corners. 
This is essentially a low-order finite element method. Although more sophisticated approximations may yield more 
more accurate strains, we are constrained by the fact that the strain is to be used with the cubic grid which is most 
convenient for the electronic structure part of the calculation. Only one strain field is needed for all sizes since it is a 
dimensionless quantity, and the strain for a different island size is obtained by a trivial rescaling. 

Because the materials lack inversion symmetry the strain produces a polarization given by P,; — CijkCTjk, leading 
to an additional electrostatic potential. For III-IV semiconductors, the only nonzero elements of the piezoelectric 
tensor are Cxyz — ^zxy — Syzx = ei4- From the polarization it is a simple matter to compute the electrostatic 
potential Vp by numerically solving Poisson's equation. Dipoles are induced at some of the edges, and on the {111} 
and {111} faces. Because the islands are primarily under compression the faces have a negative charge (Fig. 2). 
The piezoelectric effect contributes only a small amount to the potential felt by carriers in the island, and most of 
the modification of the potential is outside the island, as seen in Fig. 2. In studies of InAs islands the piezoelectric 
potential was important because the islands were assumed to have four-fold rotational symmetry which was broken 
by the piezoelectric potential The six-sided InP islands lack such a symmetry to begin with, and the piezoelectric 
potential does not make a qualitative change in the spectrum. 



III. HAMILTONIAN 

The electronic structure is calculated in the envelope approximation using atj and Vp computed above. The electron 
Hamiltonian is IfTll 



F[p — Er 



2m 



arCr — eV„ 



(3.1) 



where Ec is the local conduction-band edge in the absence of strain, Oc is the conduction-band deformation potential, 
and (T = trcTij . For the valence band we used a 6-band model given by 
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R = [l2[dl - dl) - 2*73^.9,] 

S = -V373 — dz{d.j: - idy) 
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(3.5) 
(3.6) 
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where Ey is the unstrained local valence-band edge, and ji are the Luttinger parameters. Since the Luttinger 



parameters vary spatially, we used the prescription 7^9^ 
The strain dependent coupling is given by |^] 
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(3.9) 
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where ay, b and c are the hydrostatic and two shear deformation potentials. The four band model may be obtained 
by taking A — > cxi. 



IV. MATERIAL PARAMETERS 



The values used for the various material parameters are given in table 1. These vary considerably in the accuracy 
with which they have been measured, with the deformation potentials suffering the most uncertainty. There is a paucity 
of data on GalnP, and most quantities are determined by interpolating between InP and GaP. The conduction-band 
effective mass in GalnP is simply set to the value for InP, since interpolation is complicated by the fact that GaP is 
indirect. The hydrostatic deformation potentials for GalnP, Oc and ay, are estimated by assuming ad ay is the same 
as for InP, while = a„ + Oc is taken from GalnP measurements |8||9[]. In computing Vp, the dielectric constant is 
assumed to take on the InP value throughout the system. 

One parameter that bears special mention is the unstrained valence-band offset, (i.e. the InP valence-band energy 
referenced to that of the GalnP, in the absence of strain.) The value used is based on transition-metal impurity 
spectra, and is in agreement with the value based on Au schottky barrier data |]l0|| . The idea is that transition-metal 
impurities are empirically observed to have energy levels fixed with respect to the vacuum, relatively independent of 
their host environment. Thus, by comparing band edges referenced to the impurity levels in two different materials 
one deduces the relative band offsets if the strain could be turned off. The ground state energies of Fe impurities are 
0.785 eV and 0.74eV above the valence band in InP and GalnP respectively [|ll|, so the InP valence band is 45meV 
below GalnP. 
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V. ELECTRONIC STRUCTURE 



Before solving Schrodinger's equation it is instructive to examine the band structure shown in Fig. 3. The 
conduction-band is rather ordinary, with electrons confined to the island by a barrier of about 250meV at the 
InP/GalnP interface. Strain reduces the barrier height, but since most of the band offset resides in the conduction- 
band, there is still a substantial barrier. There is also some modification of the band edge inside the island, but for 
the most part it looks like a particle in a box. The barrier height of 250meV is in good agreement with measurements 
using deep- level transient spectroscopy, which give a barrier heigth of 240meV [ [L2[ . 

The valence band has a very complex structure. Note that in the absence of strain the island would actually be a 
weak antidot since the InP unstrained valence-band offset is — 45meV; The hole confinement is due entirely to strain. 
The valence-band edge has peaks near the bottom corners of the island which extend around the base in a ring. The 
highest points in this ring lie near the {111} planes. Because of the shallowness of the hole potential, the piezoelectric 
potential has a greater impact on the holes than on the electrons. Looking along the z axis we see the valence band 
also has peaks in the GalnP immediately above and below the island. 

The confined state energies and wave functions are found by diagonalization of a finite difference version of the 
Hamiltonian using the Lanczos algorithm. This is done on the same cubic grid used to compute the strain, the only 
difference being that it is truncated to exclude unnecessary regions of barrier material. The symmetry of the island 
is not used to reduce the grid. 

The first few conduction-band wave functions are shown in Fig. 4. These states are relatively simple to understand 
if one crudely approximates the truncated pyramid as a flat box which is smallest in the z direction. The low-lying 
states all correspond to excitations in the x and y directions. One manifestation of the strain is the ordering of the 
first and second excited states which appears reversed from the naive expectation. If the island were a simple box, 
the state with a node along the [110] direction would be higher since the island is shorter in that direction. However, 
strain modifies the the potential, resulting in a different ordering of states. 

The valence-band states have a more complex structure, as seen in Fig. 5. The states fall into two categories: states 
localized near the bottom of the island, labelled and those localized above the island in the GalnP, labelled 
The ground state is Aq, which is peaked around the band-edge maxima near the {111} planes. Excited states (lower 
valence-band energy) extend around the ring at the base of the island. The first type B state, i3o, appears 15meV 
from the ground state. Although the band diagram shows a pocket below the island, there are no localized states 
there. The localized valence-band states provide a partial explanation for the observation of multiple lines seen in 
photoluminescence (PL) jl^ . An excited hole that relaxed into a B state would be unable to transfer into one of the 
A states, and would instead recombine with an electron in its ground state. Examination of the valence-band wave 
functions shows that including the split-off band is not a superfluous addition. The split off component is typically of 
order 0.2 of the largest component. 

Since island size depends on growth conditions it is interesting to examine the dependence of energy on island size, 
as shown in Fig. 6. The conduction-band spacings are of order lOmeV for a height of 15 nm, increasing to 20meV 
for a 10 nm island. The hole spacings are much smaller due to their higher effective mass, and the dependence on h 
is more complex due to the odd shaped potential. It should be noted that although the spacings re not monotonic in 
/i, the individual energies are. 

VI. TRANSITIONS 

Knowledge of transition rates are important for determining which states will be experimentally accessible. In 
order for the B states to produce an observable PL line the B ^ A transition rate must be sufficiently slow that 
recombination occurs before relaxation into an A state. A simple method of analysis is to compare the relaxation 
rates for Bq An with An An-i, where An is the yl-state closest in energy into which the Bq state could 
relax. If, for example, we assume the transition is due to emission of an acoustic phonon, the rate is given by 
Ti^f cx |(/|exp(i(f • The calculated wave functions give TBo-,An ^0~^^Ta„^a„-i for q throughout the 

Brillouin zone. This suppression is due entirely to the highly localized nature of the states, so a comparable suppression 
will be obtained for other relaxation processes. 

Since multiple PL lines are observed experimentally, it behooves us to examine optical transitions involving excited 
states. The strength and polarization dependence of the radiative transitions are found by computing the dipole 
matrix elements. These are computed by decomposing the envelope wave functions into linear combinations of the 
states {X\, {Y\, {Z\, and {S\, and using the fact that only {X\px\S) = {Y\py\S) = {Y\pz\S) is nonzero. The matrix 
elements are calculated for light in the z direction for both possible polarizations. The results are shown in Fig. 7, 
where Ix and ly are the squares of the dipole matrix elements for light polarized in the x and y directions respectively. 
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Casual examination of the wave functions might lead one to believe that the transition involving the B state would 
be suppressed since it is type II. This concern is seen to be unfounded, and the transition rate from the B state 
is comparable to that from the A state. The B transitions are highly suppressed for conduction-band states that 
are antisymmetric about the [110] direction. The A transitions are less selective, but still show strong polarization 
dependence. The most notable feature of transitions between the conduction-band ground state and the valence-band 
states An is that they are predominantly polarized in the x direction. Therefore, even if there is mixing of the A 
states, the resulting transition should still be polarized in the x direction. 



VII. CONCLUSIONS 

Comparison with PL measurements shows good agreement with the major features |l^. PL of single quantum dots 
shows energies in the range 1.62 — 1.64eV, which corresponds to the calculated energies for island heights in the range 
14 — 17 nm. The calculated PL spacing of 15meV for h = 15 nm is in agreement with the observed multiple PL lines 
which are 10 — 20meV apart. Experiments on single islands show as many as four lines, which the current model 
is unable to explain. It should be noted, however, that macroscopic PL measurements which necessarily average 
over a range of island sizes and shapes, typically show two peaks spaced 15 — 20meV apart. Therefore the dominant 
characteristic of the spectrum is consistent with the valence-band double-well potential. The additional lines appear 
to be due to some further structure. The extra structure could be due to asymmetry of the islands, which would each 
split each of the A and B states into two localized states, giving a total of four lines. While appealing, this model is 
difficult to reconcile with the observation that the islands have a high degree of symmetry ^ . 

In summary, we have shown that InP islands have a rich electronic structure, with holes confined to multiple 
pockets in and around the island. The energy levels are consistent with measurements on single dots, and the 
calculated structure can produce two PL lines 15meV apart. While the predicted doubling of PL lines provides an 
important existence proof, more work is needed to fully explain the puzzling proliferation of lines seen in experiments. 
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TABLE I. Material parameters. Unless otherwise noted, InP values are taken from Reference 11, and GalnP values are 
interpolated between InP and GaP values from Reference 11. 



Parameter 


InP 


GalnP 


nic 


0.077 mo 


0.077 mo 


71 


4.95 


5.24 


72 


1.65 


1.53 


73 


2.35 


2.21 


Eg 


1.424 eV 


1.97 eV 


A 


0.11 eV 


0.095 eV 




-6.6 eV 


-7.1 eV ^ 




-7.0 eV " 


-7.5 eV 


a„ 


0.4 eV 


0.4 eV 


b 


-2.0 eV 


-1.9 eV 


d 


-5.0 eV 


-4.75 eV 


ei4 


0.035 C/m^ " 


0.068 C/m^ ° 


EH 


12.61 


12.61 ^ 


Cxxxx 


10.22 X 10" dyne/cm^ 


12.17 X 10" dyne/cm^ 


Cxxyy 


5.76 X 10" dyne/cm^ 


6.01 X 10" dyne/cm^ 


Cxyxy 


4.6 X 10" dyne/cm^ 


5.82 X 10" dyne/cm^ 


a 


0.58687 nm 


0.56532 nm 



Value for InP used. 

Reference IBj 

Reference H. 
"^see text. 
" Reference jl|]. 

FIG. 1. Island geometry. The distances to the {011} and {111} planes measured from the center of the base of the island 
are in the ratio doii/rfin ~ 1. 

FIG. 2. Piezoelectric effect for an island with /i = 15 nm. (a) Piezoelectric charge density. The red surface is the contour 
of p = +0.001e/nm'^, blue surface indicates p = — O.OOle/nm'^. (b) Electrostatic potential due to the piezoelectric polarization, 
Vp. The red surface is the contour of = +20mV, blue surface indicates Vp — —20mV. 

FIG. 3. Band diagrams with inclusion of strain. Scale shown is for h — 15 nm, but other island sizes are obtained by simply 

rescaling. (a) Along z direction, through the center of the island, (b) Along the line x = y passing through the base of the 
island. 

FIG. 4. Conduction-band wave functions and corresponding energies for /i = 15 nm. red/blue indicates positive/negative 

FIG. 5. Valence-band wave functions for h — 15 nm. The contours are of \ipi'^ equal to 0.01 of the peak value, (a) Blue 
surface is the ground state, ^o, yellow surface is the highest state located above the island. Bo- (b) blue surface is the excited 
state An immediately above Bo in energy. 

FIG. 6. Energies as a function of island size, (a) Electron states, (b) Valence band. Only the first few A states are shown, 
so there are missing states between the lowest A state and the highest B state (-Bo)- 

FIG. 7. Magnitude and polarization dependence of band-to-band optical matrix elements for an island with h = 15 nm. 
Ix + ly is in equal arbitrary units in all three figures, (a) Transitions from first six conduction-band states C„ to valence-band 
state Aq, (b) Transitions from first six conduction-band states C„ to valence-band state Bo, (c) Transitions from conduc- 
tion-band ground state Co to valence-band states Ao...A-i4,Bo. Transitions spaced less than O.lmeV apart have been combined. 
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